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Analysis  systems  for  the  mesoscale  thermal 
variability  in  the  Greenland- Iceland- Norwegian 
Sea 

D.  Grillaki-Steiert  and  V.  Amoroso 


Abstract:  Mesoscale  thermal  structure  and  its  spatial  variability  are  very 
significant  for  acoustic  propagation  in  the  ocean.  Therefore  ocean  ther¬ 
mal  analysis  procedures  ate  important  for  the  initialisation  and  verification 
of  environmental  and  acoustical  models.  This  study  involved  the  imple¬ 
mentation  of  software,  and  the  production,  from  irregularly-spaced  data,  of 
uniformly-gridded  temperature  fields  for  use  by  numerical  acoustical  and 
environmental  models.  The  two  analyses,  considered  for  the  production 
of  temperature  fields  are:  (1)  an  objective  analysis  method  based  on  the 
Gauss-Markov  Theorem  which  makes  use  of  additional  information  such 
as  climatology  (the  Levitus  climatology  has  been  used  in  this  study),  and 
(2)  an  interpolation  method  based  on  contouring  techniques  provided  by 
the  UN  IRAS  graphics  package.  \Data  from  the  Summer  '86  cruise  in  the 
Greenland-Iceland-Norwegian  SeV  of  the  Applied^Oceanography  Group  of 
SACLANTCEN  are  used  to  demoiistraTeThe  methodSv^Both  methods  agree 
well  when  there  is  a  high  density  of  data.  When  the  data  availability  is 
limited,  there  is  a  danger  that  the  objective  analysis  method  can  give  rise 
to  spurious  gradients  in  going  from  the  regime  dominated  by  data  to  that 
dominated  by  climatology.  Whether  this  is  less  desirable  than  a  field  based 
on  extrapolation  from  a  few  observational  data  points  depends  on  the  ap¬ 
plication. 


Keywords:  acoustic  propagation  o  Gauss-Markov  theorem  o 
Greenland  Iceland- Norwegian  Sea  a'  Levitus  climatology  A  modelling  o 
oceanography  d  thermal  analysis  o  thermal  variability 
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1.  Introduction 

Oceanographic  and  acoustic  numerical  models  require  a  synoptic  (instantaneous) 
picture  of  the  ocean  thermal  field.  The  development  and  operational  use  of  range- 
dependent  acoustic  models  is  dependent  on  the  information  provided  by  such  an 
ocean  thermal  field.  However  the  simultaneous  collection  of  the  uniformly-spaced 
oceanographic  samples  required  for  the  synoptic  description  of  a  thermal  field  would 
be  too  expensive  to  realize  in  practice.  Satellite  data  could  be  used  to  provide 
a  synoptic  picture  of  the  surface  thermal  field  but  they  cannot  directly  give  any 
information  on  the  thermal  field  below  the  surface. 

In  this  work  an  objective  analysis  method  and  a  contouring  analysis  method  have 
been  used  to  transform  irregularly-spaced  data  to  uniformly-gridded  fields  for  use 
by  numerical  oceanographic  and  acoustic  models.  This  involved  two  distinct  tasks. 

(1)  The  implementation  of  software  to  adapt  the  objective  analysis  system  to  the 
Greenland-Iceland-Norwegian  Sea  (GIN  Sea)  area,  and  to  interface  it  with 
an  oceanographic  database  and  the  commercially- available  graphics  package 
UNIRAS.  Also  to  interface  the  UNIRAS  contouring  system  to  the  oceano¬ 
graphic  database. 

(2)  The  production  of  two-dimensional  high-resolution  (0.5°  x  0.5°)-gridded  ther¬ 
mal  data  fields  of  the  GIN  ’86  data  (i.e.  data  from  the  Summer  ’n6  cruise  of  the 
SACLANTCEN  Applied  Oceanography  Group  in  the  GIN  Sea)  with  the  two 
methods,  for  input  into  oceanographic  or  acoustical  range- dependent  models. 
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2.  Data 


2.1.  SAMPLED  DATA  SETS 

The  data  used  in  this  study  are  92  CTD  profiles  obtained  using  a  Neil  Brown  Instru¬ 
ment  Systems  Conductivity-Temperature- Depth  profiler.  Temperature  data  sets,  lor 
depth  levels  of  10,  50,  100,  150  and  200  m,  were  extracted  from  an  oceanograhic 
database  and  used  to  initialise,  together  with  the  Levitus  climatology,  the  objec¬ 
tive  analysis  programs,  or  were  interfaced  directly  with  the  UNIRAS  package  for  the 
production  of  two-dimensional  high- resolution  thermal  fields. 

A  data  table,  or  data  set,  is  created  after  a  database  search  that  has  imposed  some 
selection  criteria  on  the  stored  data  files. 

For  the  CIN  Sea  data  set  display  a  formatted  column  representation  was  used:  the 
first  two  columns  contain  the  longitude  and  latitude  of  the  sampled  data,  and  the 
remaining  columns  give  the  required  oceanographic  physical  parameters. 


2.2.  LEVITUS  CLIMATOLOGICAL  DATA 

The  Levitus  climatology  (1982)  is  a  synthesis  of  all  temperature  salinity  and  Oj  data 
available  at  NODC  (National  Oceanographic  Data  Centre)  analysed  in  a  consistent 
objective  manner  on  a  1°  x  1°  grid  between  the  surface  and  ocean  bottom  at  stan¬ 
dard  ocean  analysis  levels.  Annual  and  seasonal  analyses  have  been  computed.  The 
analyses  comprised  all  data  regardless  of  the  year  of  observation  (running  from  the 
late  1940’s  to  the  first  quarter  of  1978).  The  spring  Levitus  temperature  climato¬ 
logical  data  for  our  area  of  interest  have  been  used  together  with  the  CIN  ’8t>  data 
to  initialise  the  objective  analysis  programs. 
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3.  Database 

The  two  methods,  objective  analysis  and  contouring  analysis,  used  92  CTD  profiles 
(Fig.  1)  averaged  over  10  m  intervals  from  the  Summer  ’86  cruise  in  the  Greenland- 
Iceland-Norwegian  Sea.  An  oceanographic  environmental  database  has  been  de¬ 
veloped  at  SACLANTCEN  in  the  Applied  Oceanography  Group  (AOG)  for  the 
storage  of  oceanographic  and  meteorological  data  collected  by  SACLANTCEN  or 
other  oceanographic  institutes.  Before  ocean  environmental  data  are  loaded  into 
the  AOG  database  they  are  submitted  to  data-processing  by  a  system  that  checks 
acquisition  and  data  quality  (see  Fig.  2a).  Extraction  of  stored  data  sets  and  other 
operations  on  the  database  are  made  by  a  database  search  (see  Fig.  2b),  which  fol¬ 
lows  some  selection  criteria,  and  results  in  a  collection  of  files  containing  the  required 
oceanographic  data  sets. 

After  a  database  search  two  options  can  follow: 

•  Operations  on  data  files,  changing  existing  data  sets. 

•  Data  table  generation,  i.e.  generation  of  data  sets.  (This  option  has 
been  interfaced  with  the  UNIRAS  graphics  package.) 
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4.  Interpolation  metliods 


4.1.  OBJECTIVE  ANALYSIS  METHOD 

Objective  analysis  theory,  sometimes  referred  to  as  optimal  estimation,  was  first 
developed  for  application  to  meteorological  data  by  Gandin  (19(53).  The  objective 
analysis  method  used  here  is  a  variant  of  the  optimum  interpolation  method  devel¬ 
oped  by  Bretherton  Davis  and  Fandry  (1976)  for  ocean  forecasting  anti  has  been 
modified  for  the  GIN  Sea.  The  method  interpolates  irregularly-spaced  observations 
to  a  regular  grid.  It  is  a  powerful  statistical  interpolation  technique  that  uses  cal¬ 
culated  field  statistics,  i.e.  cross-correlation  between  data  points,  to  calculate  an 
estimated  uniformly' gridded  field  extending  to  locations  where  few  or  no  data  exist. 
This  estimated  field  is  such  that  it  minimises,  in  a  least-squares  sense,  the  difference 
between  itself  and  an  ensemble  of  the  observed  fields.  Given  N  observations  of  data 
Dx  at  different  locations  a  simple  linear  estimate  of  the  true  field  Fq  at  some 
other  point  x0»  which  can  be  a  grid  point,  is 

N 

Co  =  fi>  +  ]T'My.  f.),  (i) 

t=  i 

where  Fq  is  the  estimated  field  at  grid  point  r0l  f„  is  a  first-guess  tieid  at  grid  point 
rB  (i.e.  the  value  of  the  climatological  field  at  a  grid  point),  and  f,  is  a  first-guess 
field  at  points  (i.e.  the  values  of  the  climatological  field  at  the  measured  data 
positions). 

The  a,  are  coefficients  and  are  derived  using  the  following  scheme.  According  to  the 
Gauss-Markov  theorem  tiie  minimisation  of  the  expected  scpiared  error  EB  is 

El  =  «F0  -  F0)2),  (2) 

where  (.  .  .)  is  the  statistical  expectation  or  ensemble  mean. 

Considering  the  deviation  of  the  above  fields  from  the  first -guess  or  climatological 
field  when  d,  -  D,  -  f , ,  /o  =  Fq  -  f0: 

El  =  ([(Fo  -fo)  -  (F0  fo)j2> 

N 

i  =  l 

N  N  N 

=  </o>  -  25Z“,'(/o<*.->  +  X  X  a,aj(d,dj)  (3) 

1=1  1=1  J  \ 
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Minimising  Eo  with  respect  to  a,,  dE0/da,  =  0,  gives  a  set  of  N  linear  equations 
for  ait 

N 

=  (f'dj),  (4) 

«  =  1 

which  can  be  used  to  calculate  the  a,,  if  the  statistical  properties  of  the  field  (didj) 
and  (fad j)  are  known.  The  (d,d})  is  the  covariance  of  observational  data  with  itself 
and  (fodj)  the  covariance  between  the  true  field,  at  a  grid  point,  and  the  observa¬ 
tions.  These  are  functions  of  distance  and  can  be  estimated  from  the  correlation 
function  of  the  observations,  or  can  have  some  imposed  functional  form.  In  our 
study,  in  line  with  experience  gained  at  NORDA  (P.  May,  personal  comunication), 
they  have  been  given  the  form 


(didj)  =  exp  [-(a  |  xt  -  x,  I  +b  |  y{  -  |  +c  |  -  z,  |] ,  (5a) 

(fodj)  =  exp  (-(a  |  *0  -  z,  \  +b  |  y0  -  Vj  I  +c  |  z0  -  Zj  |] .  (5b) 


The  constants  a,  b  and  c  are  chosen  based  on  the  data  and  a  physical  understanding 
of  the  processes  involved.  In  our  case  we  resort  to  external  observations  and  estimate 
these  values  to  be  a  =  (50  km)-1,  b  =  (50  km)-t,  c  =  (0.5  km)-1.  They  are  the 
reciprocals  of  eddy  length  scales  in  x,  y,  and  z  direction  of  our  area  of  interest. 

Equation  (4)  is  of  the  form 

CX  =  £>,  (6) 

where  C  is  the  positive-definite  symmetric  matrix  with  the  correlations  of  all  ob¬ 
servational  pairs  (Eq.  (5a)),  D  is  the  vector  with  the  correlations  between  the  grid 
points  and  observational  data  (Eq.  (5b)),  and  X  contains  the  values  of  the  coef¬ 
ficients  which  minimise  the  error,  which  are  to  be  determined.  For  the  inversion 
of  matrix  C  we  used  a  Cholesky  factorisation  followed  by  forward  and  backward 
solution  (for  more  details,  see  Grille ki ,  1988). 

The  associated  error  field  e0,  for  the  estimated  field  Fa,  is  given  by  the  equation 

N 

e0  =  l  -  (7) 

i= i 


The  objective  analysis  method,  including  a  full  derivation  of  these  equations,  is  given 
in  greater  detail  by  Gandin  (1963),  and  Bretherton,  Davis  and  Fandry  (1976). 

The  objective  analysis  system  reads  the  observational  data  and  the  first-guess  field. 
The  observational  data  read  were  provided  by  the  92  CTD  profiles  from  the  GIN  ’86 
cruise  data,  and  the  LevitU3  climatology  was  used  as  a  first-guess  field.  An  0.5°  x  0.5° 
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grid  from  9°W  to  1°  W  longitude  and  from  60°N  to  63.5°N  was  chosen  to  generate  in¬ 
terpolated  temperature  values  at  the  grid  points  using  the  objective  analysis  method. 
The  chosen  area  has  a  high  density  of  data  and  quite  a  strong  signal  resulting  from 
the  Arctic  front  (see  Fig.  3). 


4.2.  UNIRAS  INTERPOLATION  METHOD 

The  UNIRAS  software  package  applies  three  interpolation  methods  to  data  sampled 
irregularly  in  space  in  order  to  obtain  data  on  a  regular  grid.  The  UNIRAS  interpola¬ 
tion  method  does  not  use  additional  data  sources  such  as  climatology.  The  method 
is  a  combination  of  double-linear,  quadratic  and  distance- weighting  interpolations. 

The  double-linear  interpolation  interpolates  the  four  closest  data  points  to  the  grid 
point  to  be  estimated  (see  Fig.  4).  Firstly,  the  interpolated  values  between  nearest 
neighbours  are  computed  at  lines  of  constant  *,  x  —  xo  and  y,  y  =  y0 : 

*/.  =  2d,  +  A(zrfl+i  -  zdi),  (8) 

y/.  =  Vd.  +  A(!M,+,  -  SM.),  (9) 

xl,  =  *d,  +  A(u,+1  -  xdJ,  (10) 

where  A  =  (x0  -  xd,)/(*d,+l  -  *d.)  and  i  rotates  through  1  to  4  cyclically. 

Secondly,  the  pairs  of  values  at  x  -  Zq  and  y  =  yo  are  linearly  interpolated  to  ( x0 , 
y0)  and  the  mean  value  z*“  determined: 


If  data  points  are  found  in  only  three  quadrants  (see  Fig.  5)  then  only  two  intersec¬ 
tions  with  the  x  —  Xo  and  y  =  yo  plane  can  be  calculated.  Assuming  these  t  be  z/. 
and  z/j  the  grid  node  is  assigned  the  value  z"1^: 

=  zh  +  -  yn  )•  ( 12> 


If  data  points  are  found  in  only  two  adjacent  quadrants  then  only  one  intersection 
with  the  *  =  xo  or  y  =  y0  plane  can  be  calculated.  This  value  is  assigned  to  the  grid 
point.  If  data  points  are  found  in  two  opposite  quadrants  or  only  in  one  quadrant 
then  the  nearest  data  point  is  assigned  to  the  grid  node. 

The  quadratic  interpolation  provides  another  estimate  of  the  field  but  one  with 
smoothly-changing  gradients.  Let  R\  be  the  distance  from  the  grid  point  to  the 
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closest  data  point  in  grid  cell  units,  R2  =  1.25/?i,  and  i,  and  ;2  be  the  average 
value  of  data  points  at  grid  distances  i? j  and  R2.  Then  z""  is  the  estimate  value 
at  grid  node  (z,  y),  for  see  Eq.  (11): 

!""  =  R]zt  -  R\z1/R\  - R],  (13) 

=  <t  +  0.15(zinl-  zl'dv).  (14) 

The  constants  1  25  and  0.15  are  default  factors. 

Finally  the  distance-weighting  interpolation  uses  an  inverse  distance  function.  The 
weight  of  the  grid  point  in  the  centre  is  proportional  to  the  l/dr  of  the  four  surround¬ 
ing  grid  points.  In  the  current  version  of  UNIRAS  the  parameter  r  =  3  (UNIRAS, 
1986). 

The  UNIRAS  interpolation  method  reads  the  GIN  '86  observational  data  and  gen¬ 
erates  interpolated  values  at  the  grid  points  using  each  of  the  methods  described 
above,  and  these  are  then  simply  averaged  to  produce  the  final  estimate  of  the  grid 
point  value  (for  more  details  see  UNIRAS,  1986).  The  resulting  regular  temperature 
fields  are  located  on  a  0.5°  x  0.5°  grid  from  9°W  to  1°W  longitude  and  from  60°N 
to  63.5°N  latitude. 
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5.  Results 

Following  software  implemention,  high-resolution  regularly-gridded  maps  of  the  GIN 
’86  temperature  data  were  derived  using  the  two  different  analysis  methods  described 
above.  From  the  AOG  database,  where  the  92  CTD  profiles  of  the  GIN  ’86  data 
were  stored,  the  temperature  data  sets  were  extracted  for  depth  levels  of  10  m, 
50  m,  100  m,  150  m,  200  m.  These  depth  levels  were  considered  representative  for 
constructing  the  mesoscale  thermal  field  of  the  upper  oceanic  layer  of  the  area  of 
interest.  This  area  was  from  -9°W  to  -1°W  longitude  and  60°N  to  63.5°N  latitude. 
The  measured  temperature  values  are  displayed  at  the  corresponding  CTD  locations 
as  colour  dots  in  Figs.  5a-9a.  (A  colour  scale  of  the  temperature  values  is  included 
in  each  of  the  figures.)  The  Levitus  climatological  temperature  data  interpolated  on 
a  0.5°  x  0.5°  grid  by  cubic  splines  are  shown  in  Figs.  5b-9b.  In  all  the  selected  depth 
levels  we  notice  south-north  structured  gradients  that  reveal  the  southern  edge  of 
the  Arctic  front. 

Comparing  the  GIN  ’86  data  (Figs.  5a-9a)  and  the  Levitus  climatological  data 
(Figs.  5b-9b),  we  notice  marked  differences  in  the  thermal  variability.  The  CTD 
data  points,  plotted  as  individual  dots  to  highlight  the  sampling  distribution,  reveal 
the  presence  of  a  cold  intrusion  from  the  north  ( T  <  6.5  °C  at  100  m;  see  the  blue 
dots  in  Fig.  7a)  into  the  relatively  warm  North  Atlantic  waters  ( T  >  7.0  °C  at 
100  m).  This  is  absent  in  the  climatology  (Fig.  7b).  Such  intrusions  and  eddies 
are  known  to  have  lifetimes  of  only  a  few  weeks  and  dimensions  of  50-150  km,  and 
also  suffer  from  interannual  variability,  so  that  a  climatology  representing  a  seasonal 
average  of  many  years  over  a  1°  x  1°  area  cannot  reproduce  them.  Also  at  100  m 
depth,  the  CTD  data  show  warm  Atlantic  water  (yellow  dots  in  Fig.  7a)  extending 
northwards  at  the  eastern  edge  of  the  domain,  i.e.  along  the  edge  of  the  Shetland 
Shelf,  whereas  the  climatological  data  show  this  warm  water  (T  >  8.5  °C  at  100  m; 
Fig.  7b)  failing  to  reach  the  southeastern  corner  of  the  domain. 

We  may  also  note  that  the  entire  CTD  survey  took  approximately  3  weeks  to  com¬ 
plete,  so  strictly  speaking  it  is  not  a  instantaneous  sample  of  the  thermal  field. 
Where  the  measurements  have  been  repeated  after  a  couple  of  weeks  along  the  same 
track  running  from  6°  to  2°W  and  60.5°  to  62°N  (roughly  along  the  diagonal  of  the 
plot),  significant  variations  can  already  be  noticed  (see  Figs.  5a-9a). 

The  results  of  the  data  interpolation  using  the  UNIRAS  method  are  given  in  Figs.  5c- 
9c  at  depths  of  10,  50,  100,  150  and  200  m.  The  corresponding  results  with  the 
objective  analysis  method  are  given  in  Figs.  5d-9d.  From  these  figures  we  see  that, 
in  regions  where  there  is  a  high  density  of  data  the  agreement  between  the  two 
methods  is  good.  However,  in  data-sparse  regions,  such  as  the  northwest  corner  of 
the  domain,  to  the  north  and  west  of  the  Faeroe  Islands  (~  62°N,  ~  6°W),  there 


-  8  - 


SACLANTCEN  SM-212 


are  significant  differences.  Again,  using  the  100  m  data  as  an  example  (Figs.  7c 
and  7d),  the  objective  analysis  method  produces  predominantly  zonal  contours, 
whereas  those  from  the  UNIRAS  method  are  compressed  into  a  meridional  feature  at 
the  centre  of  the  northern  edge  of  the  box.  The  objective  analysis  method,  making 
heavy  use  of  climatology  (see  Fig.  7b)  where  the  data  is  sparse,  thus  spreads  out 
the  cold  water,  T  <  6.5°,  over  the  whole  northern  quarter  of  the  region,  whereas 
the  UNIRAS  method,  relying  on  weighted  extrapolation,  describes  a  ‘tongue’  of  cold 
water  between  5°W  and  2°W  at  all  depth  levels.  The  southernmost  tip  of  the  narrow 
tongue  (defined  by  T  <  5.5  °C  at  100  m  depth)  reaches  ~  62.4°N  in  the  UNIRAS 
analysis  results,  but  reaches  only  ~  63.0°N  as  a  broad  feature,  in  the  objective 
analysis  results.  The  situation  at  100  m  depth  in  the  southeastern  corner  is  also 
markedly  dissimilar  in  the  two  representations  The  UNIRAS  method,  following  the 
CTD  data  (Fig.  7a),  makes  the  8.5  °C  contour  meander  to  the  northeast  leaving 
the  domain  at  61.5°N,  1°W,  whereas  the  objective  analysis  method,  following  the 
climatology  (Fig.  7b),  has  it  going  due  south  along  the  2°W  meridian 

When  the  data  density  is  low  the  objective  analysis  method  is  dominated  by  the  cli¬ 
matology  in  use.  The  UNIRAS  method  in  data-sparse  areas  extrapolates  a  smoother 
field  from  the  data. 
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6.  Discussion  and  conclusions 

Two  interpolation  methods  have  been  presented  in  this  report  that  produce  regularly- 
gridded  maps  of  a  field  (in  this  case,  temperature)  for  the  Greenland-Iceland- 
Norwegain  Sea  area  for  use  by  numerical  environmental  and  acoustical  models.  The 
two  methods  are  based  on  different  mathematical  principles  and  one  uses  an  addi¬ 
tional  source  of  information  (climatology)  but  their  products  agree  for  the  temper¬ 
ature  field  in  areas  of  dense  data  coverage. 

The  primary  factor  affecting  the  quality  of  both  analysis  methods  seems  to  be  the 
quantity  of  data  available.  In  practice  the  quantity  of  data  available  for  an  opera¬ 
tional  analysis  will  usually  be  quite  limited.  The  use  of  additional  information  helps 
to  overcome  this  problem,  and  the  objective  analysis  method  has  the  ability  to  utilise 
additional  independent  information,  e.g.  climatology,  to  produce  an  improved  anal¬ 
ysis.  An  important  factor  for  the  good  performance  of  the  objective  analysis  theory 
is  the  choice  of  climatology.  The  results  presented  here  also  show  that  when  the 
large-scale  structure  of  the  data  and  of  the  climatology  do  not  agree  (i.e.  when  the 
field  sampled  by  the  data  has  significant  departures  from  the  climatological  mean 
field)  the  objective  analysis  method  can  give  rise  to  anomalous  horizontal  gradients 
in  going  from  the  regime  dominated  by  data  to  that  dominated  by  climatology.  Such 
an  example  can  be  seen  in  the  south-east  corner  of  our  area  of  interest  at  100  m 
depth. 

In  oceanography  there  is  always  the  problem  of  inadequate  data  sampling.  Generally 
in  an  observational  data  set  there  will  be  large  data  gaps  or  significant  elapsed  time. 
However,  although  the  observational  data  coverage  in  this  analysis  (92  CTD  profiles 
sampled  over  a  period  of  3  weeks  covering  an  area  from  9°  to  1°  W  and  60°  to  63.5°N) 
is  denser  than  for  many  oceanographic  applications,  a  definite  comparison  between 
the  two  methods  cannot  be  made.  In  areas  of  inadequate  data  coverage  the  output 
of  the  methods  gives  rise  to  discrepancies.  The  application  will  dictate  which  type 
of  error  or  uncertainty  can  be  tolerated  (i.e.  whether  spurious  gradients  are  more 
tolerable  than  a  fictitious  extrapolated  field).  In  our  experience  the  objective  analysis 
scheme  cannot  be  applied  to  all  situations  without  exercising  some  degree  of  caution 
in  the  interpretation  of  the  output  field.  A  more  definite  study  should  be  made  in 
the  future  using  perhaps  a  synthetic  field  or  even  high-resolution  satellite  imagery. 

We  also  plan  to  compare  similar  analyses  using  different  climatologies,  such  as  the 
Generalized  Digital  Environmental  Model  (GDEM)  climatology,  and  to  develop  a 
climatological  database  from  all  AOG  cruises  and  other  available  data  in  the  GIN 
Sea  area. 
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Fig.  2a.  Data  acquisition  and  database  plan. 


Fig.  2b.  Database  plan  diagram. 


-  @2  - 


Fig.  3.  Location  of  the  Arctic  Front  over  the 
Iceland-Faroe  Ridge.  The  nid-gradient  isotherm 
envelope  at  300  m  (solid),  at  100  m  (dotted)  and 
the  mean  position  of  the  3°  isotherm  at  300  m 
(dashed). 


Fig.  4.  UNIRAS  interpolation  method  grid  scheme. 
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Fig  5  (»)  Temperature  measurements  at  10  m  depth  from  the  GIN  '86  data,  (b)  Temper¬ 
ature  field  obtained  from  the  Levitus  climatology  at  10  m  depth. 
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Fig.  5  coat'd.  Thermal  analysis  resulting  from 
objective  analysis  interpolation  using  the  GIN 
depth 
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(c)  the  UNIRAS  interpolation  and  (d)  the 
’86  data  and  Levitus  climatology  for  10  m 
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Fig.  6  cont’d.  Thermal  Analysis  resulting  from  (c)  the  UNIRAS  interpolation  and  (d)  the 
objective  analysts  interpolation  using  the  GIN  ’86  temperature  data  and  Levitus  climatol¬ 
ogy  for  50  m  depth. 
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Fig  7.  (a)  Temperature  measurements  at  100  m  depth,  from  the  GIN  '86  data  (b)  Tem¬ 
perature  held  obtained  from  the  Levitus  climatology  at  100  m  depth. 
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Fig.  7  cont'd.  Thermal  analysis  resulting  from  (c)  the  UN1RAS  interpolation  and  (d)  the 
objective  analysis  interpolation  using  the  GIN  ’86  data  and  Levitus  climatology  for  100  ni 
depth. 
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Fig.  8.  (a)  Temperature  measurements  at  150  m  depth,  from  the  GIN  '86  data,  (b)  Tem¬ 
perature  field  obtained  from  the  Levitus  climatology  at  150  m  depth. 
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Fig.  8  cont’d.  Thermal  analysis  resulting  from  (c)  the  UN  IRAS  interpolation  and  (d)  the 
objective  analysis  interpolation  using  GIN  ’86  data  and  Levitus  climatology  for  150  m 


SaCLANTCEN  SM-212 


Unimap  200— m  temp 


8.5 

8.0 


7.5 
7.0 

6.5 


6.0 

5.5 
5.0 

4.5 


4.0 


3.5 
3.0 

2.5 


GIN  SEA  OBJ.  ANA.  LEVITUS  CUM  200  MET 

63.5 


4a  *  > 

^  £  > 


63.0 

62.5 
62.0 

61.5 
61.0 

60.5 


60.0  -P 

-9  -6  -7  -6  -5  -4  -3  -2  -1 


ABOVE 

7.5 
7.0 

as 

6.0 

5.6 
5.0 

4.5 
4.0 

3.5 
3.0 

BELOW 


B.O 

0.0 

7.5 
7.0 
6  5 
6.0 

5.5 
5.0 
4  5 
4.0 
3  5 
3  0 


Fig.  9coitt’d.  Thermal  analysis  resulting  from  (c)  the  UNIRAS  interpolation  and  (d)  the 
objective  analysis  interpolation  using  GIN  ’86  data  and  Levitus  climatology  for  200  m 
depth. 
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